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ABSTRACT 

The manner in which gas accretes and orbits within circumnuclear rings has direct implications for the star formation process. In 
particular, gas may be compressed and shocked at the inflow points, resulting in bursts of star formation at these locations. Afterwards 
the gas and young stars move together through the ring. In addition, star formation may occur throughout the ring, if and when the gas 
reaches sufficient density to collapse under gravity. These two scenarios for star formation in rings are often referred to as the ‘pearls 
on a string’ and ‘popcorn’ paradigms. In this paper, we use new Herschel PACS observations, obtained as part of the KINGFISH 
Open Time Key Program, along with archival Spitzer and ground-based observations from the SINGS Legacy project, to investigate 
the heating and cooling of the interstellar medium in the nearby star-forming ring galaxy, NGC 4736. By comparing spatially resolved 
estimates of the stellar FUV flux available for heating, with the gas and dust cooling derived from the FIR continuum and line 
emission, we show that while star formation is indeed dominant at the inflow points in NGC 4736, additional star formation is needed 
to balance the gas heating and cooling throughout the ring. This additional component most likely arises from the general increase in 
gas density in the ring over its lifetime. Our data provide strong evidence, therefore, for a combination of the two paradigms for star 
formation in the ring in NGC 4736. 
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1. Introduction 

Circumnuclear star formation ri ngs are pres ent in ap proximately 
20% of nearby disk galaxies (IComeron et al.11 2010) . Early the¬ 
ories suggested that the radius at which the rings form was 
the inner Lindblad resonance (ILR) radius, which is set by an 
asymmetric component in the stellar disk such as a large scale 
bar. However, further modeling of galaxy gravitational poten¬ 
tials has shown that circumnuclear ring formation is linked more 
particularly to the pres ence of a sufficient number of *2 orbits 
dRegan & Teubenll2003l) . Such orbits are (near-)circular, in con¬ 
trast to the elongated x\ orbits atjarger radii, which make up 
the large scale bar ( Contopoulos & Papavannopoulosll980b . The 
radii of *2 orbits are close to the ILR radius. 

As gas and dust lose angular momentum, they move inward. 
At the overlap between x\ to %2 orbits, gas and dust will settle 
on the lower energy X 2 orbits as a result of collisions, forming a 
circumnuclear ring. The discrete locations where this transition 
happens are generally called ‘inflow points’. Such inflow points 
are most obvious in ring systems with large scale stellar bars, 


where the connection between the spiral shocks of the bar and 
the circumnu clear ring can often be clearly seen (for example 
NGC 6951. lYan der Laan et ai1l201 lh . However, circumnuclear 
rings also exisHn galaxies classified as non-barred (see, e.g., 
[Grouchy et al.ll201Q[ for a discussion on the difference, or lack 
thereof, between barred and non-barred systems in this context). 
Any gravitational asymmetry in the disk may suffice, as long as 
enough *2 orbits are present in the gravitational potential. 

The gas and dust can experience shocks as they transi¬ 
tion from x\ to X 2 orbits at the inflow points. However, af¬ 
ter the gas and dust settles on the X 2 orbits of the ring, the 
accreted ISM experiences very little large-scale friction. With 
that in mind, it is not too surprising that it has been found in 
multiple galaxies that most (if not all) star formation in cir¬ 
cumn u clear rings takes pl a ce at the inflow poi n ts (iRvder et al. 


2001 

; Allard et al. 2005& Boker et al. 2008; Mazzuca et al. 

2008 

; Van der Laan et al. 2013^fbj; Falcon-Barroso et al. 2014]). 


Consequently, it is thought that, owing to the absence of large 
scale friction or temporal variations in gravitational potential, 
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stars and gas primarily move through the circumnuclear ring 
together. The stellar age gradient regularly seen stretching out 
from each inflow point is often referred to as ‘pearls-on-a- 
string’, first coine d bv lBoker et al.l (120081) . An alternative theory 
(Elmegreen 1994) suggests that the inflow points play no spe¬ 
cial role in star formation, and that a build up of gas will trigger 
a burst of star formation after a critical limit has been reached 
anywhere in a ring. This is often referred to as the ‘popcorn’ 
theory of star formation in rings. 

Recent star formation is often tracked by observing hydrogen 
recombination lines such as H a or Bry. The equivalent width of 
the Baimer li nes is, for example, dependent on stellar age from 
~4 Myr (e.g. LeithererHeckman 119951: Bresolin & Kennicutt 
1 19971: iGonzalez Delgado et al.l ll999h . Alternatively, emission 
lines from other atoms and molecules stimulated by stellar 
radiation, such as warm H 2 , Hel, and [Fell], in the optical 
and NIR can be used. With a combination of these lines it 
is also possible to observe potential age gradients, as demon¬ 
strated for the circumnuclear rin g in NGC 613 dBoker et al.l 
20081: iFalcon-Barroso et al.l [20141) . However, strong extinction 
can make observations in the optical and NIR range difficult, 
especially at the inflow points where dust and gas build up, po¬ 
tentially limiting the effectiveness of these methods. 

When there are no large-scale mechanical stresses, stel¬ 
lar radiation is most likely to be the prevalent form of heat¬ 
ing for the neutral gas in the ring (excluding contributions of 
CR or X-ray ionization). Thus, an alternative way of track¬ 
ing recent star formation is to study the heating and cooling 
of the neutral interstellar medium (ISM) via strong emission 
lines of the neutral gas in the FIR. In the presence of star for¬ 
mation, the heatin g of the neutral ISM i s dominated by photo¬ 
electric heating dBakes & Tielensl 11994 IWeingartner & PraTnel 
120011: iHelou et al.ll200ll) . Dust grains and PAHs eject electrons 
as the result of the absorption of FUV (6eV < hv < 13.6eV) 
photons from young O and B type stars. These electrons in turn 
interact with the gas, heating atoms and molecules. The cooling 
times for dust and gas are much shorter than the Myr astronom¬ 
ical timescales of interest here. Thus, gas heating and cooling, 
via photo-electric heating, may be considered in balance at all 
times, and the observed cooling flux can be used as a proxy for 
the heating. 

To test whether the heating of the neutral ISM in the ring is 
indeed dominated by stellar radiative heating from stars formed 
at the inflow points, we investigate the case of the circumnuclear 
ring in NGC 4736, a nearby (D=4.66Mpc, i ~ 35°, PA~ 300°) 
SAab spiral galaxy. The ring in NGC 4736 has a radius of 
40" (0.9 kpc). Based on a previously observed HI velocity field 
([Walter et al l 120081) and assuming that the spiral arms are con¬ 
cave, the north side of the disk is the near side, and the south 
the far side, with a clockwise rotation. We have combined avail¬ 
able data to form a large multiwavelength dataset, tracing the 
gas and the recent star formation, including Herschel PACS 
data from the KINGFISH (Key Insights on Nearby Galaxies: 
A Far-Infrared Survey with Herschel) Open Time Key Program 
dKennicutt et al J120111) . The PACS observations enable us to 
probe the main neutral gas cooling lines of [CII] at 158yum and 
[OI] at 63pm, as well as the FIR dust continuum. We also in¬ 
clude data obtained with the Spitzer Space Telescope as part of 
the SINGS (Spi tzer Infrared Nearby Galaxies Survey) Fegacy 
program dKennicutt et al.120031) . to probe the mid-infrared wave¬ 
length range, which includes emission from large dust grains, 
PAH’s, stars, and molecular and atomic gas. 

In §2 we present the data. In §3 the atomic and molecular 
gas mass surface densities, and subsequent dust mass surface 



Fig. 1 . Spatial coverage of the spectral data sets presented in this 
study, overlaid on the GALEX FUV image of NGC 4736. The cir¬ 
cumnuclear ring at a radius of r~40" (0.9 kpc) stands out clearly 
in the GALEX image. The different regions mapped spectro¬ 
scopically, and indicated here in color, are: Spitzer/IRS long-low 
(14-38 pm) spectra (dark blue contours), Herschel/PACS [CII] 
and [OI] emission line maps (green contours), Herschel/PACS 
[Nil] 122yum, [NII]205/im emission line maps (yellow contours), 
and Spitzer/IRS short-low (5-14/mi) spectra (purple contours). In 
this image, north is up and east is to the left, and the total size is 
approximately 7.3 arc minutes on a side. 


density, are computed. In §4 a census is made of the star for¬ 
mation rate density throughout the ring. In §5 the contribution 
of the neutral ISM to the observed flux is determined. We then 
compare the dust and gas cooling with the stellar heating in §6. 
The implications are discussed in §7, with a summary in §8. 


2. Data sets 

2.1. Herschel PACS 

Both spectra and images were obtained for NGC 4736 with the 
Herschel PACS instruments as part of the KINGFISH Open 
Time Key Program. The PACS spectroscopic instrument is an 
IFU (integral field unit) with 5x5, 974 pixels dPoglitsch et al.l 
?()!(» . This gives an instantaneous Field-of-View (FoV) of 
47" x47". The PACS beam FWHM is 572 in the wavelength 
range 60-85/mi, and 12" at 130-210yum. For the spectra sev¬ 
eral targeted pointings were made for the various lines. The 
[NII]122yum and [NII]205yum lines were observed with 2 adja¬ 
cent pointings, one focussed on the nucleus and one on the cir¬ 
cumnuclear ring (see yellow contours in Fig. GO). The [CII]158/mi 
and [OI]63yum lines were observed in a strip 200" long (green 
contours), with an orientation equal to previous Spitzer IRS data 
(blue contours), and one adjacent pointing on the circumnuclear 
ring (equal to the [Nil] off-nucleus pointing). The FIR contin¬ 
uum was probed at 70yum, 100/mi, and 160/mi with the PACS 
photometer over the full extent of the galaxy. 

The PACS observations were reduced consistently for the 
complete KINGFISH sample using the Scanamorphus package 
(Roussel 20131). The KINGFISH observa tions are described in 
iDale et al.l (120121) and lCroxall et al.l (120121) . In all cases, the im¬ 
ages were rescaled to a pixel scale of 4". The spectra zeroth- 
moment maps were convolved to the PACS 160/im resolution of 
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Fig. 2. Distribution of the Total InfraRed (TIR) continuum emission (left), [CII] 158 pm line emission (middle), and [01] 63 p m line 
emission (right) in the circumnuclear ring of NGC 4736. The 12" apertures used to extract fluxes are indicated in the TIR and [CII] 
panels. The spatial scale is given in the bottom right corner, and the PACS 160/rm beam of -12" is indicated in the top left corner 
of each panel. 


12" , using a pixel scale of 4" , since the PACS 160 beam is es¬ 
sentially the same as the [CII] beam. 

The [CII] and [01] intensity maps, as well as the total-IR 
(TIR) continuum map are shown in Fig. [2] The TIR continuum 
is constructed from a combination of four individual Spitzer and 
Herschel images, at 8 pm, 2Apm, ID pm, and 160 pm, following 
the equation in lDraine & Lil (2 0071 , 

TIR = (0.95(v5' y ) 8 + 1.15(v 5 v >24 + <vS y ) 7 o + <vS y ) 16 o) • (1) 

Even though the FoV of the spectra clips the northern and south¬ 
ern part of the circumnuclear ring, the [CII] and [01] emission 
is clearly concentrated there. The TIR continuum also has ele¬ 
vated emission at those locations, but is mostly dominated by 
the nucleus. 

Integrated fluxes are extracted by covering the ring with 3x3 
pixel (12"xl2") apertures. In this manner we measured the 
[CII] and [01] emission line and continuum fluxes at 13 indi¬ 
vidual points, seven on the east side, and six on the west side of 
the ring. For the TIR continuum an additional six apertures are 
positioned on the ring outside the spectroscopic FoV (see Fig. [2]). 
This east (west) side corresponds to the blue (red) shifted side of 
the galaxy. The positions of, and fluxes in, these 12" apertures 
are given in Table [T] 

2.2. Spitzer 

As part of the SINGS survey NGC 4736 was observed with 
Spitzer -IRS using the ‘short-low’ and ‘long-low’ spectral mod¬ 
ules in spectral mapping mode. The circumnuclear ring is in¬ 
completely covered by the ‘long-low’ footprint (see again Fig.Q] 
blue contour), and partially by the three off-nuclear ‘short-low’ 
pointings (purple contours). The ‘short-low’ (SL) module has a 
pixel scale of 1785 and a resolving power of 60-120, decreasing 
with increasing wavelength. The observations have a FoV of 376 
by 27" and cover the wavelength range 5-14/rm. The ‘long-low’ 
(LL) module has a pixel scale of 571, a resolving power of ap¬ 
proximately 60-120, again decreasing with wavelength, a FoV of 
1076 by 84" , and covers the wavelength range 14-28yum. The re¬ 
duced and fully calibrated data cubes are available on the SINGS 
Spitzer and Ancillary Data websiteB 


1 http://irsa.ipac.caltech.edu/data/SPITZER/SINGS/ 


For this work we extracted the [Sill] emission line flux at 
18.7 pm and 33.5 pm by fitting the lines and the underlying warm 
dust continuum with PAHFIT ([Smith et al .1120071) on the LL data 
cube alone. The [Sill] emission lines are not blended, and the 
dust emission is a smoothly varying composition of black-body 
curves of different temperatures. Consequently, the separation of 
line and continuum emission at these wavelengths with PAHFIT 
is straight-forward. A check at the positions with both SL and 
LL data showed that the exclusion of the SL data cubes had no 
detrimental effect on the fitting of the long-wavelength dust con¬ 
tinuum. Since the IRS data are not combined with the other data 
sets - they will only be used in Sect. [5]- and have a limited FoV, 
no effort was made to convolve the cubes to the PACS 160 reso¬ 
lution. 

Spitzer-IRAC imaging at 8 pm and 24 pm are also presented, 
and were convolved to the PACS 160 resolution using kernels 
from lAniano et al] (1201 ll) . 

2.3. Ancillary data 

The THINGS surve)0 (IWalter et alJl200S) observed the HI dis- 
tribution of NGC 4736 with the VLA between March 2003 and 
March 2004 in the BCD configurations. The HI intensity maps 
(natural and robust weighting) are publicly available from the 
THINGS website. The beam size of the observations is, de¬ 
pending on the chosen weighting, 10722 x 9707 (natural), or 
5796x5755 (robust). __ 

The BIMA SONG survey dHelfer et al .1 2 00 3) observed the 
CO(l-O) distribution in NGC 4736 in October 1997. The total 
FoV of these data are -190" (4.3 kpc). The data have a reso¬ 
lution of 6" and are short spacing corrected. The resulting data 
cube and moment maps are publicly available from NED. 

NGC 4736 was observed with GALEX in both the NUV and 
FUV bands between April and June, 2004. The observations 
are publicly available through the MAST GALEX public access 
site|j. The FoV is 1.2°, with a FWHM angular resolution of 470 
(NUV) and 576 (FUV). 

An H a map was obtained from the NAS A/IPAC Extraglactic 
Database (NED), but was originally obseiwed, reduced, and con¬ 
tinuum and [Nil] subtracted bv lKnapen et al.1 ( 2004 ). It was ob- 


2 http://www.mpia-hd.mpg. de/THINGS/Overview.html 

3 http://galex.stsci.edu/GR6/ 
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Table 1. [CII], [01], and TIR surface brightnesses in defined apertures 


# 

pos X 
[pix] 

pos y 
[pix] 

Angle 

[deg] 

[CII] 

W/m 2 /sr 

[CII ] U nc 
rel. 

[OI] 

W/m 2 /sr 

[OIW 

rel. 

TIR 

W/m 2 /sr 

TIR WWC 

rel. 

1 

149 

144 

15 

- 

- 

- 

- 

1.29e-4 

0.09 

2 

146 

143 

29 

1.12e-6 

0.01 

2.43e-7 

0.12 

2.64e-4 

0.07 

3 

143 

145 

49 

1.13e-6 

0.01 

3.33e-7 

0.06 

2.84e-4 

0.04 

4 

141 

148 

68 

1.62e-6 

0.01 

4.91e-7 

0.04 

4.54e-4 

0.04 

5 

141 

151 

85 

1.75e-6 

0.01 

5.40e-7 

0.03 

4.88e-4 

0.03 

6 

142 

154 

102 

1.60e-6 

0.01 

4.24e-7 

0.05 

3.98e-4 

0.03 

7 

143 

156 

117 

l.lle-6 

0.01 

2.72e-7 

0.09 

2.94e-4 

0.03 

8 

145 

158 

135 

8.87e-7 

0.02 

2.00e-7 

0.23 

2.56e-4 

0.04 

9 

148 

159 

157 

- 

- 

- 

- 

2.70e-4 

0.03 

10 

151 

159 

180 

- 

- 

- 

- 

2.56e-4 

0.03 

11 

155 

161 

205 

8.50e-7 

0.02 

1.31e-7 

0.19 

2.20e-4 

0.03 

12 

157 

161 

213 

7.83e-7 

0.02 

1.69e-7 

0.11 

2.12e-4 

0.03 

13 

160 

160 

229 

1.05e-6 

0.01 

2.94e-7 

0.06 

2.66e-4 

0.03 

14 

161 

157 

243 

1.34e-6 

0.01 

3.70e-7 

0.05 

3.28e-4 

0.03 

15 

161 

154 

259 

1.44e-6 

0.01 

3.82e-7 

0.05 

3.54e-4 

0.03 

16 

161 

151 

275 

1.40e-6 

0.01 

4.34e-7 

0.08 

3.38e-4 

0.03 

17 

159 

148 

297 

- 

- 

- 

- 

2.88e-4 

0.05 

18 

156 

146 

320 

- 

- 

- 

- 

3.42e-4 

0.10 

19 

154 

145 

337 

- 

- 

- 

- 

3.14e-4 

0.12 

20 

152 

144 

353 

- 

- 

- 

- 

2.56e-4 

0.12 


Notes: [CII] and [01] emission line fluxes, TIR continuum flux, and their relative uncertainties in 3x3 pixel (12" xl2") apertures. Central 
positions (in TIR map pixels) of the apertures are given in the 2nd and 3rd columns. The order of the apertures (angles) here is clockwise 
(subsequently following the direction of gas motion), starting from due south. For reference, pixel (151,152) corresponds to the galaxy center 
(RA 12h50m53.18s, Dec 41d07ml2.0s). 
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Fig. 3. Left and middle: Mass surface density maps of atomic (HI) and molecular (CO(l-O)) gas, at the PACS 160yum (~12") reso¬ 
lution. Right: Total gas mass surface density map, based on HI and CO maps, with a 26% correction for helium. The radius of the 
circumnuclear ring and the locations of the inflow points are indicated with a dashed circle and stars, respectively. Clockwise gas 
flow is indicated with the white arrows. 


tained with the Isaac Newton 2.5m telescope. The angular reso¬ 
lution is 174x174 and the pixel scale is 073. 

3. ISM distribution in the ring 

3.1. Atomic and molecular gas 

Both the HI and CO(l-O) intensity distributions, convolved to 
the PACS 160/mi resolution, are shown in the first two panels 
of Fig. [3] The HI distribution is dominated by two arcs which 
are co-spatial with the circumnuclear star forming ring (whose 
radius is approximated with the dashed ellipse, representing an 
i ~ 35° inclined ring). The nucleus is devoid of atomic gas. By 
comparison, the CO(l-O) distribution shows emission in a more 
spiral arm-like pattern. (This is even more apparent at native res¬ 
olution.) The spiral arms only partly overlap with the circum¬ 


nuclear ring. Unlike the HI, the CO(l-O) intensity distribution 
further shows strong emission from the nucleus. 

The ‘total’ gas mass density is computed via these two trac¬ 
ers, by assuming a CQ-conversion factor , Xco, of 7.4 xlO 19 
[cm -2 (K km/s) _1 ] (ISandstrom et al.ll2013l) . This value is lower 
than the Galactic X C o value of 2.0 x 10 20 [c m~ 2 (K km/s)" 1 ] 
most o ften assumed for nearby galaxies, but ISandstrom et al.1 
( 2013 ) have shown that Xco drops to this value towards the 
galactic nuclear region of this galaxy (averaged at a resolution of 
3775). Due to the lower Xco factor, the molecular gas composes 
a smaller fraction of the total gas mass. The total gas mass map 
is still dominated by molecular gas, as can be expected in the 
molecule rich centers of galaxies, but atomic gas provides about 
1/3 of the gas mass in the ring. A 26% correction for helium and 
metals has been applied to the total gas map. 
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Ring inflow points are defined as the position where the gas 
spiral arm(s) cross onto the ring. In a bi-symmetric potential, e.g. 
with a stellar bar, the two points should be more or less sym¬ 
metric about the nucleus. An indication of their position in this 
galaxy is given in Fig. [3] The total gas mass map shows two, 
possibly three, regions of high mass surface density at the radius 
of the ring, in the east, the northwest, and at lower significance 
to the southwest. The first two are near the inflow points. 

3.2. Dust 

By comparing our gas mass map with a map of the dust mass, 
we can determine whether or not the dust and gas in this region 
are well mixed, which is a necessary condition for our upcoming 
analysis and comparison to models of the star formation in the 
ring. The dust mass distribution in NGC 4736 is derived using 
via the iDraine &TJ ( 2007 ) model, including Spitzer, Herschel 
PACS0 and SPIRE observations of the source, to constrain the 
IR dust emission (lAniano et al.l 120121) . We find that our mass 
distribution of gas shows a one-on-one morphology correspon¬ 
dence, confirming that dust and gas are qualitatively well mixed. 
Equally, when we compare our TIR continuum map (Fig. [2]) with 
the dust luminosity found by the model a one-on-one morphol¬ 
ogy correspondence can be observed. Additionally we find that 
our total gas mass map also matches quantitatively, i.e., with rea¬ 
sonable conversion values for the dust-to-gas ratio (DGR) we 
recover the values in the dust map. 

4. Star formation distribution in the ring 

The star formation rate in the circumnuclear ring can be derived 
from either the Her, NUV, or FUV data sets. To compensate for 
the stellar light absorbed by dust, all three star formation rate 
tracers are combined with the Spitzer 24yum image to obtain 
dust-corrected values. 

The star form ation rate is computed based on the SFR cal¬ 
ibrations given in iMurphv etakl (1201 ll) . lHao et al.l (1201 ll) . and 
Ken nicutt & Evansl (120121) . The resulting SFR surface density 
(SFRD) maps are shown in Fig.® Differences between the three 
maps reflect the uncertainties and physical/temporal variations 
among the SFR tracers. The region for which the SFRD is com¬ 
puted is, in all cases, limited at small radii by a circular mask at 
20" to exclude emission from the nucleus. At the outer edge, a 
limit has been placed at the minimum of either 70" radius, or 
the position where the flux in the tracers falls to a 2 cr level. The 
maps shown in Fig. [4] are at the resolution of the 24yum observa¬ 
tions, which is -5" , and have a pixel scale of 175. SFRD maps 
at the PACS 160 /i m resolution were also computed. 

In all three panels, the regions of highest star formation rates 
are found in the southeast and northwest quadrants of the ring. 
These positions are near/at the ring inflow points, as expected. 
At the other two quadrants of the ring, the SFRD is lower (by 
factors of three to six). This could indicate a secondary star 
formation event of lower-level intensity, or the continued emis¬ 
sion from an ageing burst. The SFRD contrast around the ring 
is strongest in the Hcr+24/im map. The SFRD valuesfound for 
Ha+24//m were checked a gainst previous results bv lWong et al.l 
(l2004l) . IWong et al.l (120041) found a SFRD for the circumnuclear 
ring, based on H a fluxes alone, of 0.04Mo/yr/kpc 2 . We match 
this value with our H a observations. 

Comparing with the gas maps shown in Fig.® the SFR dis¬ 
tribution is most similar to that of the atomic gas. This is surpris- 

4 the same observations as used in this work 


ing, given the fact that stars form from molecular not atomic gas. 
We speculate here that the young stars in the ring of NGC 4736 
are energetic enough to dissociate large quantities of molecu¬ 
lar gas in their immediate vicinity (without ionizing that gas as 
well), and will investigate this matter in an upcoming paper (Van 
der Laan, in prep.) 

The star formation regions near the inflow points each cover 
about a fifth of the ring. The orbital time for a star at this ra¬ 
dius (r=0.9kpc) is about 25Myr, given a 220km/s rotational 
velocity (value derived from the BIMA SONG first moment 
map). Therefore the duration of the starburst cannot be more 
than 5Myr. This is also a reasonable time span to explain 
the difference between the Hcr+24yum, and NUV+24/rm and 
FUV+24yum maps. When we ‘subtract’ the H a distribution from 
either the NUV or FUV distribution (not shown), what remains 
is a more uniform distribution of flux, expected if the UV emis¬ 
sion drops more slowly with time (-5.5 Myr for H a, - 11 Myr for 
NUV/FUV, as based on Starburst99 stellar burst synthesis mod¬ 
els). The UV tracers are still tracing the previous ‘event’ when 
they have moved through half the ring. 

With these observations it is not possible to directly distin¬ 
guish an age-gradient of star formation in the hot spots seen at 
the inflow points. Inspection of the observations shows that the 
hot spots are strongest in the 24 jum observations and nearly un¬ 
seen in the FUV and H a, implying they are (strongly) extincted 
regions. However, an equivalent width map of the H a emission, 
given that a suitable continuum image is found, could potentially 
show age-related variations between the hot spots. 

Averaged over the ring, the SFRD values correspond to a star 
formation rate density of -0.065 M©/yr/kpc 2 . Within individual 
hotspots, the SFRD reaches more than 0.25M©/yr/kpc 2 . High 
SFRD are a common characteristic of circumnuclear star form¬ 
ing r ings (iKormendv & Kennicuttl 120041: iKennicutt & Evansl 
120121) . Integrating over the full area of the ring leads to SFR val¬ 
ues of 0.16, 0.27, and 0.25 M©/yr, for Ha+24yum, NUV+24yum, 
and FUV+24yum, respectively. These values are lower than 
previously derived S FRs for the whole gala xy of 0.38M©/yr 
( Calzetti et al.l 20101) . 0.43M©/yr dLeel [2006). and 0.70M©/yr 
(Ski bba et al.ll2011 ). but show the major role the circumnuclear 
ring plays in the overall star formation in NGC 4736. 


5. Radiation field in the ring 

In order to investigate the effect the young stars have on the lo¬ 
cal ISM, it is important to measure the emission from the dense 
neutral regions (the PDRs) close to the forming stars. 

As a first step, it is important to estimate, and then remove, 
the contribution of [CII] from the diffuse, ionized ISM. Neutral 
Carbon has an ionization potential (11.3eV) which is slightly 
lower than that of Hydrogen (13.6eV). Consequently, it will 
be present in both the neutral and ionized ISM. The contribu¬ 
tion of [CII] from the diffuse ionized ISM can be determined 
from an observed [CII]/[NII]205 jum ratio, since N + emission ar¬ 
rises exclusively from the diffuse ionized medium dBennett et al.l 
1 994 , and the critical density of [CII] and [Nil] 205/rni are 
nearly the same. The observed [CII]/[NII] ratio is nearly con¬ 
stant (-4) in the ionized ISM over a large range of electron 
density, n e . The difference between the predicted and measured 
[CII]/[NII]205/im ratio can then be used to estimate the ionized 
and neutral fractions of [CII] emission. 

Since the [NII]205/im line is not significantly detected 
in the NGC 4736 ring, we follow an analysis similar 
to the one discussed in ICroxall et al.l ( 2012 ) and use the 
[CII]/[NII]122//m ratio to determine the ionized [CII] emission. 
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Fig. 4. Star formation rate density maps, based on Ho'+24 y um (left), NUV+24yum (middle), and FUV+24/mi (right), at the Spitzer 
24/rm (~5") resolution. The Hcr+24/im derived SFRD values in this image have been multiplied by a factor two for visualization 
purposes. The galactic nucleus and larger radii have been masked. The radius of the circumnuclear ring and the locations of the 
inflow points are indicated with a dashed circle and stars, respectively. The Spitzer 24yum beam of ~5" is indicated in the top left 
corner of each subpanel. 
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Fig. 5. Plot of [CII]/[OI] vs. ([CII]+[OI])/TIR to determine 
the ambient UV radiation field incident on the PDR regions. 
The [CII] emission has been corrected for a 10% contribu¬ 
tion from the diffuse ionized ISM. The underlying grid is the 
[Kaufman et al. ( 199 9) PDR model for different densities and 
stellar radiation field strengths. The inset sub-panel is a zoom- 
in of the region of interest. Data points in the east/west side of 
the ring are indicated in blue/red. 
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Fig. 6. Same as Fig. [5} except with a further assumption of 20% 
TIR, and 60% [CII] contribution from PDRs illuminated by a 
high radiation field intensity, with the rest originating from more 
diffuse neutral ISM components. 


The [CII]/[NII]122/im ratio is dependent on electron density in 
the range 1-100 cm -3 . While the [Sill] 18.7/33.5yumratio is avail¬ 
able from our IRS maps, the ratio [Sill] 18.7/33.5 is sensitive to 
much higher electron densities in the range 10 2-5 cm -3 , and the 
[Sill] ratio (0.56) is consistent with gas in the low density limit. 

The [Nil] 122/205 ratio is sensitive to electron densities of 
0.1-100 cm -3 . Since the [Nil] 205/rm line is not detected in 
the ring, the observed [NII]122yum signal-to-noise ratio in the 
ring (5 - 10) can be used as a proxy for the upper limit of 
the [Nil] 122/205 ratio. From lBeirao et alJ (120121) we know that 
[Nil] 122yum/[NII]205yum increases with increasing electron den¬ 
sity, with the ratio reaching 10 at n e = 100 cm -3 . Thus, in the 
ring we are dealing with an n e lower than 100 cm -3 , consistent 
with our Spitzer/IRS measurements. At an n e of 100 cm -3 the 
[CII]/[NII]122yum ratio from the diffuse ionized component of 
the ISM would be 1.24 (this is an upper limit). All additional 
[CII] emission should be from the neutral component. In the re¬ 
gion of the ring for which we have [NII]122yum emission maps, 
we find [CII]/[NII]122/im ratios of 8-16, with the value increas¬ 
ing away from the highest intensity regions. These values corre¬ 
spond to 85 to 92% contributions of the neutral phase to the [CII] 
emission. At lower electron density the [CII] fraction in the neu¬ 
tral phase rises. We therefore assume that 90% of the measured 
[CII] emission in the ring is from the neutral phase. 

In Fig.[5] the (diffuse ionized emission corrected) [CII]/[OI] 
and ([CII]+[OI1)/TIR ratios are overlaid on a PDR model from 
[Kaufman etal. ( 1999 , hereafter K99). As described in Sect. 12.11 
the [CII], [OI], and TIR fluxes were extracted in square apertures 
of 12" by 12" size, covering the circumnuclear ring. No differ¬ 
ence is seen between the two (Fig.0 red/blue points) sides of the 
ring. Based on these [CII]/[OI] and ([CII]+[OI])/TIR ratios the 
radiation field Go, a dimensionless variable in multiples of the 
radiation field in the Solar neighborhood, in the PDR regions of 
the circumnuclear ring is found to have a uniform value of ~200, 
and the gas density, n H , in the PDR regions to be ~ 10 3 cm -1 . 

We can compare our Go derived above from the FIR emis¬ 
sion lines with the average <U> obtained from modeling the 
dust emission. <U> is a measure of the int erste llar radiation field 
in units of the solar neighborhood dDraine & Lill2007l) . For NGC 
4736, <U> is roughly 10 (Aniano et al., priv. comm.), about a 
factor of 20 lower than the Go inferred above. This large mis¬ 
match between the radiation field intensity derived from the gas 
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Fig. 7. Integrated FUV flux between 6-13.6 eV in the case of an 
instantaneous starburst, rescaled to 1M©. Within 4Myr the en¬ 
ergy has halved, and within 11 Myr the energy has dropped by 
an order of magnitude. (Based on Starburst99 v6.0.3. models) 


and dust has been noted in other studies as well dCroxall et al.l 
20121) . and can be understood by noting that the PDR fraction, 
the fraction of the dust illuminated by the high intensity radiation 
field, can be much less than unity. Aniano et al. find a f pdr (the 
fraction of dust mass illuminated by Go greater than -110) in 
the ring of order 20%. If we apply this fraction to the TIR, and 
assume only 60% of the neutral [CII] emission is from dense 
PDRs illuminated by the high radiation field intensity, we can 
reconcile the two estimates (see Fig. 

6. Heating-cooling balance 

In the ‘pearls-on-a-string’ model, the heating in the ring as a 
function of azimuth is greatest close to the inflow points, fol¬ 
lowed by a slow decay as the stars age and the ring rotates. To 
see if the observed cooling in the ring follows this pattern, we 
compare the dust continuum flux and the far-IR line emission, 
against the expected available stellar heating. 


6.1. Energy input into the ISM 

By making several reasonable assumptions, the available en¬ 
ergy for dust and gas heating from the newly formed stars in 
the ring can be derived from the measured SFRD. An instanta¬ 
neous starburst was modeled using Starburst99 dLeitherer et al.l 
119941201 (Th . assuming a standard Kroupa IMF and solar metal- 
licity. The output synthetic line spectrum was integrated between 
912 A and 1120 A (6-13.6eV), to obtain the photo-electric heat¬ 
ing energy as a function of time. The resulting available heating 
power per M© is shown in Fig. 0 As can be seen, the power 
in the 6eV to 13.6 eV range of the integrated stellar population 
spectrum is strongly dependent on time. The flux plateaus in the 
first 3 Myr, but within 11 Myr, it has dropped by an order of mag¬ 
nitude. 

Instead of a single instantaneous starburst, physically the star 
formation will occur over a longer time interval. Like a Ferris 
wheel, the ring can be thought of as having discrete volume el¬ 
ements (‘cabins’), scooping up gas over a fixed time interval as 
it enters onto the ring. In this model, each volume element will 
have a star formation event lasting a fixed fraction of the rotation 
time, with the mean age of the event depending on the distance 
from the inflow point. By tracing back the observed star forma¬ 
tion under this premise, we can recreate the heating flux through¬ 


out the ring. As a reminder: the average SFRD measured in the 
ring was 0.065 M©/yr/kpc 2 . To be able to compare the available 
heating to the observed TIR emission later on, we multiply by an 
average 24yum/UV flux ratio of 0.33, based on our observations, 
as a proxy of how much FUV flux is heating the dust, rather 
than immediately escaping the galaxy. This ratio is close to the 
measured dust/stellar flux ratio of 0.25 bv lSkibba et al.l ( 201 ll) . 


6.2. Comparison to dust 

A SFR history (for one full 25 Myr orbit) is postulated where 
there are two starbursts of 5 Myr duration, 12.5 Myrs apart, with 
a burst SFR intensity of 0.01M o /yr (thus a total of 5xl0 4 M© 
stellar mass per starburst is created). The heating curve is mod¬ 
eled by convolving the instantaneous starburst, shown in Fig. 
0 with this proposed (box-)function of SFR history to obtain 
the instantaneous starburst heating as a function of time/position 
throughout the ring (left panel Fig. 0 grey saw-tooth curve). 
This heating curve is compared with the TIR continuum sam¬ 
pled throughout the ring (colored points). We may assume an 
instantaneous response from the dust to the heating, since the 
timescales here are of order -Myr, while dust responds to heat¬ 
ing on order 10 4 yr, regardless of dust geometry (B. Draine, priv. 
comm.). In Fig.0the correspondence between the measured TIR 
continuum and the modeled heating is very good. However, it is 
clear that the heating power falls too steeply around the ring, 
and that the second peak is significantly lower than the first. A 
much better match (black curve) can be made by shortening both 
initial bursts, lowering of the second burst to 0.008 M©/yr, and 
adding two short bursts of star formation at ~8.75-10Myr and 
20-21.25Myr. The inclusion of the latter is a divergence from 
the ‘pearls-on-a-string’ scenario. 


6.3. Comparison to PDR gas 

To compare the response of the PDR gas to the modeled heat¬ 
ing we need an observational measure of the photoelectric heat¬ 
ing efficiency. This is the ratio of photon energy going into gas 
heating to the total available energy. This ratio, typically taken 
as ([CII] + [OI])/TIR, is equal to the ratio of gas cooling line 
flux to TIR continuum flux. In Fig. 0 this gas to dust ratio is 
plotted against, the ratio (vS v >7o/(vS v )ioo, which is used as a 
proxy for the dust temperature. From Fig. 0 it can be seen that 
the photoelectric heating efficiency shows very little spread over 
the circumnuclear ring, with values of 0.004-0.006, about 0.5%. 

In the right panel of Fig.0 this fixed 0.5% scaling has been 
applied to the [CII]+[OI] flux in each region. It is immediately 
clear that the spread in the emission line flux is due to the FUV 
heating variation. The two cooling lines again follow the mod¬ 
eled FUV heating quite closely, although we lack the data points 
away from the inflow points to confirm the secondary star form¬ 
ing events, due to limited Herschel spectral coverage. 

The photoelectric heating efficiencies are comparable to the 
values s een for other kpc-siz ed star forming region s, as for ex¬ 
ample in lCroxall et al] (1201 2l) or lParkin et al.l(l2013l) . These stud¬ 
ies find values of -0.5% at lower dust temperatures, with a sig¬ 
nificant d rop of heating efficiency at hi gher dust temperat ures. 
However, iLebouteiller et al] (120121) and [Parkin et alJ (12014 do 
not find a reduction in the gas heating efficiency at higher dust 
temperatures in a star forming region of the LMC, and the disk 
of Centaurus A, respectively. The NGC 4736 ring does, however, 
have a higher photoelectric heating efficiency, at only marginally 
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Fig. 8. Left: Available FUV heating flux (black saw-tooth line) as a function of azimuth, derived from an proposed star formation 
history. The proposed star formation history is indicated at the bottom of the panel (black box function). The grey lines (underlying 
the black) are an alternative example of a star formation history and FUV heating flux. The observed FIR (dust) flux is shown in 
blue/red (east/west side of the ring). Details in the text. Right: Similar, except now for the [CII]+[OI] (gas) emission lines. The FoV 
of the Herschel spectra does not cover the full ring, thus for several apertures the gas information is lacking. 



<uF u 70>/<uF^00> 

Fig. 9. Gas heating efficiency proxy, ([CII] + [OI])/TIR, as a 
function of the dust temperature (blue (east]_ and red ( wes t) 
dots). For comparison, the results from IBeirao et~aD ( 20121 ) 
in the circumnuclear star forming ring of NGC 1097 (‘x’), 
and ICroxall et al.l (12012l) in sub-kpc regions in NGC 1097 and 
NGC 4559 (dashed line) are indicated. The [CII] has been cor¬ 
rected for its 10% diffuse ionized component. 


higher dust temperature, than starburst ring in NGC 1097 ring 
(iBeirao et al.llll2012l) . 


7. Implications for ring star formation theory 

We started this investigation to study if, following the ‘pearls- 
on-a-string’ paradigm for circumnuclear star forming rings, all 
star formation takes place at the inflow points of the ring. In 
this scenario, as the system rotates, the subsequent star forming 
events should create a clear azimuthal age gradient dBoker et al l 
120081 ) . 

Our simple modeling of the PDR dust and gas heating as 
function of azimuth requires secondary star formation events to 
be consistent with the observed cooling flux. While the exact 
strength, timing, and duration of these secondary events in the 
NGC 4736 ring is model dependent, it is clear that their absence 
leads to a too large contrast in heating/cooling flux around the 
ring. Since the FUV heating power available from young stars 


diminishes by an order of magnitude in 11 Myr, it seems unlikely 
that the drop be compensated for by stellar ’after-glow’. 

Are these ‘secondary’ star forming events truly random, tak¬ 
ing place from left-over gas from the ‘primary’ event? Or, are 
both events part of the ‘popcorn’ scenario (lElmegreenl 119941) . 
where star formation occurs anywhere in the ring? Given the 
relative strength of the events at both inflow points we con¬ 
sider this unlikely. This seems consistent with the fact that 
a number of galaxies show strong star formation events at 
the inflow points of their circumnuclear rings (iRvder et al] 


200 i ll Allard et al 120051: iBoker et al .1120081: iMazzucaet al. 2008; 

Van der Laan et al.ll20 13 alibi). 


A combination of ‘pearls-on-a-string’ and ‘popcorn’ seems 
like a reasonable scenario to explain the star formation in the ring 
of NGC 4736. We know star formation is not a 100% efficient 
process. Star formation efficiencies (SFE) are usually estimated 
to be of order a few percent. Although there may be an indica¬ 
tion that the SFE in the center of NGC 4736 is actually higher 
(10%, K. Sandstrom, priv. comm.), even then a significant por¬ 
tion of gas remains. Therefore, we can ask if the remaining gas 
has enough chance to cool and is dense enough to collapse again. 

The total gas mass surface density in the ring remains 
~50-80M o /pc 2 away from the inflow points (Fig. 0). This is 
higher than the ~8M 0 /pc 2 star formation threshold noticed in 
iBigiel et all (2 0081) . Additionally, while the temporal distance 
between the inflow points and these modeled secondary bursts 
is short (a few Myr), it is comparable to free fall times of giant 
molecular clouds. 


Can we expect this hybrid star formation scenario in all cir¬ 
cumnuclear rings? The life times of circumnuclear rings are of 
order 1-2 Gyr. Inflow rate s of gas onto the ring of 2 M n /y r are 
a reasonable upper limit (I Van der Laan et al.ll201 ll [2013a), al¬ 
though values of 0.1-0.25 M©/yr are more likely. So, even at ex¬ 
treme SFEs and low gas inflow rates, the gas mass in a ring will 
see an increase of about 10 8 M 0 (0.1 M 0 /yrxl Gyr) over the life 
time of the ring. Given the observed range of radii and widths 
found for circumnuclear rings (40-1000 pc radii, ~r/2 widths, 
IComeron et al. ( 2010 )), this translates easily to gas mass surface 
densities similar to those found in NGC 4736 (10-10 2 M o /pc 2 ). 
Thus, the prevalence of ‘popcorn’ star formation events should 
increase with the lifetime of the circumnuclear ring, even if star- 
bursts, ‘pearls-on-a-string’, continue to be at the inflow points. 
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8. Summary 

New Herschel FIR observations of the circumnuclear ring in 
NGC 4736 have been combined with archival observations of 
the star formation, and atomic and molecular gas distributions. 
This combination of data sets is used to characterize the mode 
of star formation in the ring. Due to the circular orbits of the 
ring itself, gas and dust will only experience large-scale shocks 
and compression at the inflow points, the locations where ma¬ 
terial transitions onto the ring. Star formation events, therefore, 
should be uniquely found close to the inflow points. We have in¬ 
vestigated the heating and cooling of the interstellar medium in 
the ring of NGC 4736, under this assumption. The stellar FUV 
flux available for heating was derived from measurements of the 
star formation rate, geometry-based estimates of time scales, and 
standard assumptions on the initial mass function (IMF). The 
FIR continuum and gas emission line fluxes were measured from 
Herschel PACS and archival Spitzer observations, and used to 
estimate the heating and cooling of the dust and gas around the 
ring. 

We find that, while heating from a significant burst of star 
formation at the inflow points can explain the gas and dust cool¬ 
ing there, an additional star formation component is required at 
other positions in the ring. This additional component is smaller 
than the star formation events at the inflow points, and most 
likely arises from the general increase of gas density in the ring 
over its lifetime. Therefore, in the NGC 4736 ring, it is a combi¬ 
nation of the ‘popcorn’ and ‘pearls-on-a-string 4 models that best 
explains the data. 
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